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We determine rapidity asymmetry in the production of charged pions, protons and anti-protons for 
large transverse momentum [pr) for d+Au collisions at ^/s^m = 200 GeV. The identified hadrons are 
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measured in the rapidity regions \y\ < 0.5 and 0.5 < \y\ < 1.0 for the pt range 2.5 < pt < 10 GeV/c. 
We observe significant rapidity asymmetry for charged pion and proton+anti-proton production in 
both rapidity regions. The asymmetry is larger for 0.5 < \y\ < 1.0 than for \y\ < 0.5 and is almost 
independent of particle type. The measurements are compared to various model predictions employ- 
ing multiple scattering, energy loss, nuclear shadowing, saturation effects, and recombination, and 
also to a phenomenological parton model. We find that asymmetries are sensitive to model param- 
eters and show model-preference. The rapidity dependence of tt~ /tv^ and p/p ratios in peripheral 
d-l-Au and forward neutron-tagged events are used to study the contributions of valence quarks and 
gluons to particle production at high pT- The results are compared to calculations based on NLO 
pQCD and other measurements of quark fragmentation functions. 

PACS numbers: 25.75.-q,25.75.Dw,13.85.-t 



I. INTRODUCTION 

The mechanisms for particle production in d-f Au 
collisions at RHIC may be different at forward 
and backward rapidities. The partons from the 
deuteron-side (forward rapidity) are expected to un- 
dergo multiple scattering while traversing the gold 
nucleus. Those on the gold-side (backward rapid- 
ity) are likely to be affected by the properties of the 
nucleus. A comparative study of particle production 
at forward and backward rapidity can be carried out 
using a ratio called the rapidity asymmetry (lAsym), 
which is defined as, 

V f ^ ^b(pt) 

where Yp and Ib are forward and backward parti- 
cle yields, respectively. YAsym may provide unique 
information to help determine the relative contribu- 
tions of various physics processes to particle produc- 
tion, such as multiple scattering, nuclear shadowing, 
recombination of thermal partons, and parton satu- 
ration. 

Recently, models incorporating different physics 
effects have described the nuclear modification fac- 
tor for d-l-Au coUisions (i?dAu)- Models including 
shadowing effects or nuclear modifications to the 
nucleon parton distributions reproduce reasonably 
well i?dAu for inclusive charged hadrons I"]. Those 
based on transverse momentum broadening (Cronin 
effect 0), dynamical shadowing, and energy loss in 
cold nuclear matter 3, also give i?dAu predictions 
for inclusive charged hadrons, consistent with ex- 
perimental data. Models based on the color glass 
condensate (CGC) approach reproduce the px de- 
pendence of inclusive charged hadron i?dAu at both 
mid- and forward-rapidity jj]. These models also 
qualitatively describe the pseudorapidity asymmetry 
for inclusive charged hadrons in d+Au collisions Q • 

Another approach based on hadronization by re- 
combination of thermal partons at lower px has been 
quite successful in describing the observed i?dAu for 
charged hadrons at RHIC u\. This approach em- 
phasizes the hadronization portion of the final state 



interaction. Although it takes into account the hard 
scattering in pQCD, the fragmentation is replaced 
by recombination of soft and shower partons in the 
intermediate pt region.^ Also, a phenomenological 
approach, called EPOS 7], based on a parton model, 
has described the d-|-Au collision data at RHIC. In 
this model the nuclear effects are included through 
elastic and inelastic parton ladder splitting. 

It is of interest to see how these models com- 
pare to data for rapidity asymmetry of identified 
hadrons from rf-l-Au collisions. More precisely, iden- 
tified hadron YAsym, a more differential quantity, 
may allow some determination of the relative con- 
tribution of the physical processes discussed above. 
Strong particle type (baryon and meson) dependence 
of the nuclear modification factor and azimuthal 
anisotropy at intermediate px (2 < pt < 6 GeV/c) 
has been observed in Au-h Au collisions at RHIC &] . 
The present study will investigate if such particle 
type (baryon and meson) dependence is observed in 
YAsym for d+An collisions. 

In addition to providing insight into different par- 
ticle production mechanisms at forward and back- 
ward rapidity for c?-|-Au collisions, the measurements 
presented here may be used to study the presence 
of possible effects of valence quarks and isospin on 
particle production. At high px and rapidities away 
from midrapidity, the role of valence quarks becomes 
increasingly dominant. Such studies are even more 
interesting for n-tag events (events where the neu- 
tron in the deuteron does not interact with the gold 
nucleus). A comparative study between p4-Au (n- 
tag) and d+Au data is of interest. For n-tag events 
at forward rapidity and high px, the two valence u 
quarks in the proton of the deuteron should lead to 
more production of 7r+ (ud) compared to 7r~ (du). 
For backward rapidities, if the flavor distribution in 
sea quarks is uniform and the incoming gold nucleus 
has no asymmetry in u and d quarks, one expects the 
ratio 7r^/7r+ ~ 1. This difference between forward 
and backward rapidity may be more pronounced for 
p/p. Study of particle ratios as a function of rapidity 
at high Px in peripheral and n-tag events for d+Au 
collisions may provide some information on the fla- 
vor dependence of particle production. These ratios 
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are in principle sensitive to the fragmentation func- 
tion ratios of u-quarks to tt" and tt"*" , to the ratio 
of (u, d)-quarks fragmenting to protons |0, and to 
the fractional contributions of quarks and gluons to 
hadrons at the given momentum. 

In this paper, we present the first results for the 
rapidity asymmetry of charged pion, proton and 
anti-proton production at high for d-f Au colli- 
sions at ^/sN^= 200 GeV measured by the STAR 
experiment [l^ at RHIC. A similar study for inclu- 
sive charged hadrons has been reported in Ref. 0. 
The asymmetry is studied as a function of pt for 
different collision centralities in the two rapidity re- 
gions \y\ < 0.5 and 0.5 < \y\ < 1.0. In section II 
we discuss the detectors used in the analysis, trig- 
ger and centrality selection, particle identification at 
high Pt, and the systematic errors. In section III we 
discuss the rapidity, px, species, and centrality de- 
pendence of YAsym- In scction IV, the YAsym results 
are compared to calculations from various models 
discussed earlier. In section V, we present the ra- 
pidity dependence of the nuclear modification factor 
for TT"*" -I- 7T~ and p+p . In section VI, we study the 
anti-particle to particle ratios as a function of rapid- 
ity at high Pt in n-tag and peripheral d-|-Au events 
in order to investigate the flavor dependence of par- 
ticle production. Section VII completes this work 
with the summary of our findings. 



II. EXPERIMENT AND ANALYSIS 

A. Detectors 

For the present analysis we use data recorded by 
the Time Projection Chamber (TPC) ^ in the 
STAR experiment at RHIC. The TPC is STAR's 
primary tracking device. It is 4.2 m long and 4 m 
in diameter. The sensitive volume of the TPC con- 
tains PIO gas (10% methane, 90% argon) regulated 
at 2 mbar above atmospheric pressure. The TPC 
data are used to determine particle trajectories, mo- 
menta, and particle-type through ionization energy 
loss [dE/dx). Its acceptance covers ±1.8 units of 
pseudorapidity (77) and the full azimuthal angle. 
Charged particle momenta are determined from the 
TPC data for the d-|-Au run in the year 2003 in 
which star's solenoidal magnet field was set to 
0.5 T. Two Zero Degree Calorimeters (ZDCs) ^ 
situated along both sides of the beam axis, about 
18 m from the nominal collision point (center of 
TPC), were used for triggering. The collision cen- 
trality is obtained from the charged hadron multi- 
plicity measured by STAR's Forward Time Projec- 
tion Chambers (FTPCs) [31 ■ The details of the 
design and other characteristics of the detectors can 
be found in Ref. ^ . The details of the trigger con- 



TABLE L Centrality selection, number of participat- 
ing nucleons, and number of binary collisions for d-|-Au 
collisions at a/snn ~ 200 GeV. 



% cross section N^^^-^ 


(A^part) 


(iVbin) 


0-20 


> 17 


15.67 ± 1.07 


15.1 ± 1.15 


20-40 


11-17 


11.16 ± 1.25 


10.6 ± 1.38 


40-100 


< 11 


5.14 ± 0.47 


4.2 ± 0.51 


0-100 


> 


8.31 ± 0.34 


7.5 ± 0.38 



dition, collision centrality selection, and method of 
high Pt particle identification are described below. 



B. Trigger conditions 

The ZDC in the Au beam direction, which is 
assigned negative pseudorapidity (77), was used as 
the trigger detector for obtaining the minimum bias 
data. The minimum bias trigger required at least 
one beam-rapidity neutron in the ZDC. The trig- 
ger efficiency was found to be 95±3% of the d-|-Au 
hadronic cross section crf^adr- Trigger backgrounds 
were determined using data recorded for beam cross- 
ings without collisions. For the n-tag events, the 
ZDC in the deuteron beam direction was used. Such 
events were required to have at least one beam ra- 
pidity neutron in the ZDC. The cross section for 
such a process was measured to be (19.2±1.3)% of 
(Jhadr- The vertex was reconstructed for 93±1% of 
triggered minimum bias events. A total of 11.7 mil- 
lion minimum bias d-|-Au events and 2.0 million n- 
tag events having a vertex within ± 30 cm of the 
nominal interaction point along the beam direction 
were analyzed. Two rapidity regions were used: 
\y\ < 0.5 and 0.5 < \y\ < 1.0, and the pt range was 
2.5 < Pt < 10 GeV/c. ThepT spectra were corrected 
for trigger and vertex-finding inefhciencies. Further 
details of trigger conditions for the minimum bias 
data can be found in Ref. . 

C. Collision centrality 

Uncorrected charged track multiplicity (A^^"^^^^) 
measured within -3.8 < rj < -2.8 by the FTPC was 
used to determine the collision centrality for d-l-Au 
collisions. The centrality selection criteria is given 
in Table along with the average number of bi- 
nary collisions (A^bin) and the number of participat- 
ing nucleons {Npart) estimated using a Monte Carlo 
Glauber calculation incorporating the Hulthen 
wave function of the deuteron p^ . In this model 
cr^;^^"^=2.21±0.09 b, and N^n for n-tag events is 
2.9±0.2. This model provides reasonable agreement 
with the measured charged track multiplicity distri- 
bution in the FTPC and the single neutron cross 
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4.0GeV/c < p.,. < 4.5GeV/c~| 




-10 



■4 -2 2 4 

number of o(7i)±5 



FIG. 1: (Color online) dE/dx distribution normalized by 
pion dE/dx at 4.0 < px < 4.5 GeV/c and 0.5 < |r;| < 1.0, 
shifted by ±5 for positive and negative charged parti- 
cles, respectively. The distributions are for minimum 
bias d+Au collisions. The pion, proton, and anti-proton 
peak positions are indicated by arrows. 



are shifted by ±5 for positive and negative charged 
particles, respectively. The ncrj distribution is a nor- 
mal Gaussian distribution with an ideal calibration. 
The six Gaussians are for tt^, and p(p). The 
na^ — ncrj and na^^^^ — ncrj values are estimated 
by studying the ncrj — ncr^, ntrj — ncrP, and ncrj - 

na^ distributions. The widths of the six Gaussians 
are taken to be the same. The Gaussian distribution 
used to extract the pion yield and the pion, proton 
and anti-proton peak positions are also shown in the 
figure. 

The proton yield is obtained by integrating the 
entries {Y) in the low part of the dE/dx distribution, 
about 2.5a away from the pion dE/dx peak. The 
integration limits were varied to check the stability 
of the results. The kaon contamination is estimated 
via either of the equations given below. The raw 
proton yield is 

V^{Y~P{h~^))/{a-P) 



or, 



we briefly describe the identification procedure (see 

< 10 GeV/c, there is a 



section measured by the ZDC on the deuteron side. 
Further details of centrality tagging in d-f Au colli- 
sions can be found in Ref . [I3 . 



D. Particle identification at high px 

Particle identification at high transverse momenta 
(pT > 2.5 GeV/c) is done by exploiting the rel- 
ativistic rise of the ionization energy loss. Here 

jriefly de 
Ref. [Hill). For 2.5 < pr 
difference of about 10-20% between the pion dE/dx 
and the dE/dx for kaons and protons, due to the rel- 
ativistic rise of the ionization energy loss for pions. 
This results in a few sigma (l-3g) separation. The 
dE/dx resolution is 8% ITsf . 

Pions are the dominant component of the hadron 
yield for d+Au collisions at RHIC. The prominent 
peak in the dE/dx distribution is used to deter- 
mine the pion yield in this range. To extract 
the pion yield in a given px bin, we performed a 
six Gaussian fit to the normalized dE/dx distribu- 
tions of positive and negative hadrons simultane- 
ously. The normalized dE/dx in general is defined 



as na^ = log{{dE / dx)Y / B x) / cr x , where X,Y can 
be , 77=*= , K"^ or p{p) . Bx is the expected mean 
dE/dx of particle X, and <tx is the dE/dx reso- 
lution of the TPC. Fig. ^ shows a typical dE/dx 
distribution normalized to pion dE/dx (referred to 
as the ncr^ distribution) for charged hadrons with 
4.0 < pt < 4.5 GeV/c and 0.5 < |r/| < 1.0. For clar- 
ity of presentation, the na^ distributions in Fig. ^ 



p={Y- PKl)/a, 

where a and (3 are the proton and kaon efficiencies 
from the integration described above, derived from 
the dE/dx calibration, resolution, and the Bichsel 
function [TsL l20l | . In the first case the kaon contam- 
ination is estimated through the yields of the inclu- 
sive hadrons {h) and pions; in the second case using 
Kq measurements [13 (only available for \y\ < 0.5 
up to pt < 5 GeV/c). The typical values of a for a 
dE/dx cut slightly away from the proton peak posi- 
tion is 0.4. The (3 values decrease from 0.2 to 0.08 
with pt in the range 2.5 < pt < 10 GeV/c. At high 
Pt, the yields of other stable particles (i.e., electrons 
and deuterons) are at least two orders of magnitude 
smaller than those of pions and are negligible for our 
studies. The two results are consistent in the region 
where STAR iiTg measurements are available. Since 
the energy loss of particles in the TPC is almost in- 
dependent of charge sign, the dependence of 
on na-jr is due to different particle composition and 
the dE/dx separation between pion, kaon and pro- 
ton [1»|. This provided a consistency check for the 
yields. 

The dE/dx resolution is better for longer tracks, 
shorter drift distance, stronger magnetic field, 
smaller multiplicity and lower beam luminosity. Due 
to longer tracks and shorter drift distances for parti- 
cles produced at higher y, the dE/dx resolution gets 
better. Thus, the separations between pions and 
kaons or (anti-)protons were larger for 0.5 < \y\ < 1.0 
than for \y\ < 0.5 and particle identification is 
easier at larger pT . 
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FIG. 2: (Color online) High transverse momentum spectra (pT > 2.5 GeV/c) of charged pions, proton, and anti- 
proton for the rapidity regions \y\ < 0.5 (solid symbols) and 0.5 < \y\ < 1.0 (open symbols) for d+Au collisions and 
various event centrality classes at ^snn = 200 GeV. 



TABLE II: Correction factors for identified hadron spec- 
tra at high pt (> 2.5 GeV/c) for minimum bias d-l-Au 
collisions. 



Type 


% 


Trigger efficiency 


95 ± 3 


Vertex efficiency 


93 ± 1 


Track reconstruction efficiency 


~ 90 ± 8 


{\y\ < 0.5) 




Track reconstruction efficiency 


~ 82 ± 8 


(0.5 < \y\ < 1.0) 




Background contamination 


~ 5 ± 1 



E. Correction factors 

The various correction factors for the identified 
hadron spectra are hsted in Table ^ The trig- 
ger and vertex efficiencies were discussed previously. 
The identified hadron track reconstruction efficiency 
was estimated by embedding Monte Carlo parti- 
cles into the real data and then following the full 
reconstruction procedure. It was observed to be 
independent of pt for px > 2.5 GeV/c for both 
rapidity regions. The reconstruction efficiency for 
Pt > 2.5 GeV/c for charged pions and protons are 
~ 92% and ~ 90%, respectively, in the rapidity re- 
gion |y| < 0.5. For 0.5 < \y\ < 1.0, the reconstruc- 
tion efficiency for charged pions and protons is ~ 
82% and 84%, respectively. The background con- 
tamination in the pion spectra for px > 2.5 GeV/c, 
primarily from Kq weak decay, is ~ 5%. No strong 
centrality dependence was observed in the correction 



factors. The charged pion, proton and anti-proton 
spectra are corrected for efficiency and background 
effects. The inclusive proton and anti-proton spec- 
tra are presented without hyperon feed down cor- 
rections 10, 1^ . Preliminary study shows that the 
ratio of A to inclusive p in the rapidity range \y\ < 0.5 
decreases from 0.7 to 0.3 with increase in from 
2.5 GeV/c to 5.5 GeV/c. 

F. Systematic errors 

The total systematic uncertainties associated with 
the pion yields are estimated to be ^ 15%, and those 
for proton and anti-proton yields are estimated to be 
^ 22%. They are of similar order for both the ra- 
pidity regions, and the average values for minimum 
bias collisions are given in Table UlTI 

The sources of systematic error on the high 
yield arise owing to: (a) uncertainty in modeling the 
detector response in the Monte Carlo simulations, 

(b) momentum resolution (increases with px) [13, 

(c) difference in the yields for different TPC sectors 
as a result of spatial distortion effects, (d) uncer- 
tainty in determining the pion and proton dE/dx 
peak positions, (e) uncertainty in estimating the 
kaon contamination to proton yields (increases from 
7% at px = 2.5 GeV/c to 15% at px = 10 GeV/c), 
and (f) uncertainty due to protons from hyperon de- 
cay that are reconstructed as primordial protons at 
a slightly higher px than their true value, with a 
worse momentum resolution (increases from 2% at 
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TABLE III: Systematic errors for identified fiadron minimum bias yields at higli pt (> 2.5 GeV/c) for d+Au collisions. 



Sources of uncertainty % Error 

Modeling detector response 8 

Momentum resolution 4 (at pT ~ 7 GeV/c) 
Spatial distortion 8 
dE/dx pion peak position 8 
dE/dx proton peak position 8 

Kaon contamination to proton yield 12 (at pT ~ 7 GeV/c) 

Protons from hyperon decay 7 (at pT ~ 7 GeV/c) 
Normalization (trigger and luminosity) 10 



PT = 2.5 GeV/c to 10% at pr = 10 GeV/c). There 
is an additional 10% |l5j normalization uncertainty 
due to trigger and luminosity uncertainties. These 
systematic errors are not shown in Figure |21 

As this work focuses mainly on ratios such as 
YAsym, most of the systematic errors cancel. The re- 
sultant systematic error on lAsym is about 5%. The 
errors shown for figures with ratios are statistical 
and systematic errors added in quadrature. 

Figure 12 shows the measured invariant yields of 
charged pions, protons, and anti-protons for the pt 
range 2.5 < px < 10 GeV/c in the rapidity re- 
gions, \y\ < 0.5 (solid symbols) and 0.5 < |y| < 1.0 
(open symbols) for minimum bias and various 
collision centrality classes for d-|-Au collisions at 
^snn = 200 GeV. The pt spectra are corrected 
for the trigger, vertex and reconstruction eSiciencies 
and the background effects listed in Table HII 



III. RAPIDITY ASYMMETRY 

In this section we discuss the y, Pt, species and 
centrality dependence of lAsym- 



loss in cold nuclear matter. 

(b) The lAsym values are consistent with unity for 
both rapidity regions at high px (> 5.5 GeV/c), sug- 
gesting absence of nuclear effects on particle produc- 
tion in d+Au collisions for this px range. 

(c) lAsym for charged pions is greater than unity 
and decreases monotonically with increasing px for 
2.5 < Px < 5 GeV/c. Although YAsym forp-|-p is also 
greater than unity, the trend seems to be towards 
a constant value in this px range. These features 
are opposite to predictions from models based on 
incoherent initial multiple partonic scattering and 
independent fragmentation p2j . Such models pre- 
dict that YAsym is less than unity at intermediate px 
and approximately unity for larger px |H . 

(d) For \y\ < 0.5, ^Asym for P+P is slightly larger 
than it is for charged pions for 2.5 < px < 4 GeV/c. 
For 0.5 < \y\ < 1.0 no strong particle type depen- 
dence is observed for YAsym- This is in contrast to 
the observed baryon-meson differences for the same 
Px range for Au+Au collisions, which were described 
by recombination-based models p^ . 



A. Rapidity, transverse momentum, and 
species dependence 

Fi gure 121 shows the high px dependence of YXsym 
for 7r++7r^ and p+p for the rapidity regions \y\ < 0.5 
and 0.5 < \y\ < 1.0 for minimum bias events. The 
backward rapidity is considered as the gold-side and 
corresponds to the negative rapidity region. The for- 
ward rapidity is the deuteron-side and corresponds 
to positive rapidity. For comparison the pseudora- 
pidity asymmetry for charged hadrons is shown 
also in the figure. The following observations are 
made: 

(a) YAsym is found to be larger for 0.5 < |?;| < 1.0 
than for \y\ < 0.5 for all the hadrons with 
2.5 < Px < 5 GeV/c. This may indicate the presence 
of some rapidity dependence of nuclear effects such 
as parton saturation, nuclear shadowing, or energy 



B. Centrality dependence 

In Fig. 01 we show the centrality dependence of 
YXsym at high px for 7r+-|-7r~ and p-l-p for the 
two rapidity regions \y\ < 0.5 (left panels) and 
0.5 < \y\ < 1.0 (right panels). The data are shown 
only for 2.5 < px < 5.5 GeV/c. The YAsym val- 
ues approach unity for the centrality classes studied 
for Px > 5.5 GeV/c in both the rapidity regions. 
For \y\ < 0.5, a prominent centrality dependence of 
YAsym is not observed. For 0.5 < \y\ < 1.0, YAsym 
is larger for central (0-20%) compared to peripheral 
(40-100%) events for 2.5 < px < 4 GeV/c. The 
indication of a centrality dependence in the YAsym 
at 0.5 < \y\ < 1.0 is consistent with predictions 
from saturation models |^. However in such mod- 
els the centrality dependence is much stronger than 
observed in the present data |5|| . 
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FIG. 4: (Color online) Centrality dependence of 
high transverse momentum rapidity asymmetry factor 
(yAsym) for Tv'^ +7v~ and p+p at \y\ < 0.5 (left panels) 
and 0.5 < \y\ < 1.0 (right panels) for d+Au collisions at 
y/s^ = 200 GeV. 



FIG. 3: (Color online) High transverse momentum ra- 
pidity asymmetry factor (VAsym) for 7r^-4-7r~ and p+p 
for \y\ < 0.5 and 0.5 < \y\ < 1.0 for minimum bias d+Au 
collisions at ^/swN = 200 GeV. For comparison the in- 
clusive charged hadron results from STAR Q are also 
shown by the curves. 



IV. MODEL COMPARISON 

In this section we compare the measured high 
Pt identified hadron lAsym in the rapidity regions 
\y\ < 0.5 and 0.5 < \y\ < 1.0 for minimum bias d+Au 
colhsions at ^snn = 200 GeV with predictions from 
various models (Figs. |SHlll- 



A. Comparison to the nuclear shadowing 
model 

First we compare the high pT charged pion YAsym 
in both rapidity regions with model predictions that 
incorporate only nuclear shadowing 1] . In these cal- 
culations two parameterizations of nuclear shadow- 
ing, covering the extremes of gluon shadowing at 
low X, are taken. The parametrization by Eskola et 
al. jM] is referred to as EKS98. The other, FGS, is 
from Frankfurt, Guzey, and Strikman (FGSO, 
the original parametrization, along with FGSH and 
FGSL for high and low gluon shadowing). The 
calculations use MRST leading order (LO) parton 
distribution functions 26]. The fragmentation of 
produced partons into charged pions uses the LO 
Kniehl-Kramer-Potter (KKP) fragmentation func- 



tions |23 obtained from a fit to e++e^ data. In 
EKS98 the valence quark shadowing is identical for 
u and d quarks at the minimum momentum scale 
of the hard interaction. In FGS the EKS98 valence 
quark shadowing ratios are used as input, along with 
Gribov theory and hard diffraction. The charged 
hadron -RdAu was reasonably well-described by such 
a model using the FGS parametrization T| . 

Our charged pion data (Fig. |SJ| indicate that 
nuclear shadowing as implemented in the models 
discussed, cannot explain the measured lAsym for 
2.5 < PT < 5 GeV/c for both \y\ < 0.5 and 
0.5 < \y\ < 1.0. The differences between data and 
model increase with increasing rapidity. At larger 
Pt, the data values approach unity, indicating an 
absence of nuclear effects. The effect on lAsym of 
using a different parametrization of nuclear shadow- 
ing at high pt is found to be negligible for \y\ < 0.5. 
However, some differences are observed in FGS for 
0.5 < \y\ < 1.0. 

^sym from the nuclear shadowing model covering 
the extremes of gluon shadowing at low x is not con- 
sistent with the measured values. The comparison 
therefore provides an idea of the maximum contribu- 
tion to YAsym from only nuclear shadowing in d-f-Au 
collisions at i/snn = 200. 

B. Comparison to the multiple 
scattering+shadowing+energy loss model 

Next we compare the high pT charged pion lAsym 
to a model that includes only coherent multiple scat- 
tering, which leads to transverse momentum broad- 
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FIG. 6: (Color online) High transverse momentum ra- 
pidity asymmetry factor (FAsym) for 7r^-|-7r~ at \y\ < 0.5 
and 0.5 < \y\ < 1.0 for minimum bias d-f An collisions 
at -^snn = 200 GeV compared to models incorporating 
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nuclear matter 01 . See text for more details. 



clear matter. The calculation, labeled as Vitev-I, has 
a slightly larger effective energy loss in cold nuclear 
matter than the one labeled Vitev-II. For \y\ < 0.5, 
both the Vitev-I and Vitev-II results are in reason- 
able agreement with the data within errors. For 
0.5 < \y\ < 1.0, the Vitev-I result slightly overpre- 
dicts the measured Ya. 



sym ■ 



The model calculations 
beyond > 3 GeV/c are independent of pt, while 
the measured yAsym tends to decrease with p^. 

For the Vitev models the rapidity dependence of 
i^Asym seems to be sensitive to effective energy loss. 
The decrease in lAsym with p^ in the model is re- 
stricted to Pt < 2.5 GeV/c. It will be impor- 
tant to have predictions from this model for proton 
and anti-protons to investigate the possible particle 
species dependence of multiple scattering in d-|-Au 
collisions. 



ening (Cronin effect), and to calculations with the 
addition of power corrections (dynamical shadow- 
ing) and energy loss in cold nuclear matter Q. 
In this model a systematic calculation of the co- 
herent multiple parton scattering in p+A collisions 
is carried out in terms of the perturbative QCD 
factorization approach. It also incorporates initial 
state parton energy loss in the perturbative calcu- 
lations. We observe (Fig. O that in both rapid- 
ity regions model expectations from the Cronin ef- 
fect, for 2.5 < Pt < 5 GeV/c, are in qualitative 
disagreement with the data. This indicates that 
multiple scattering is not the source of the observed 
asymmetry. Fig. shows also a comparison of the 
charged pion data with results of calculations which 
incorporate multiple scattering, dynamical shadow- 
ing, and a varying degree of energy loss in cold nu- 



C. Comparison to the EPOS model 

In Fig. [7| the measured lAsym for 7r++7r^ and 
p+p for 0.5 < \y\ < 1.0 are compared to the results 
from the EPOS model. In the EPOS model elas- 
tic and inelastic parton ladder splitting are the key 
processes. A parton ladder refers to the dynamical 
process of parton-parton scattering with successive 
emission of partons. The emission process can be an 
initial state, space-like cascade, or final state, time- 
like cascade. The elastic splitting in this model can 
be related to screening and saturation, while inelas- 
tic splitting is related to the hadronization process. 
This phenomenological model has been very success- 
ful in describing the inclusive charged hadron d-l-Au 
data 0. 
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The EPOS model predictions (v. 1.14) are con- 



sistent with the measured 
charged pions and p+p- 



Yas 



values for both 



weak species dependence and almost similar de- 
pendence of YAsym IS predicted. 



D. Comparison to the recombination model 

The recombination model reproduces some of the 
observed features of RHIC data 0, It suc- 

cessfully describes the Cronin effect for d-j-Au data 
without any need for fcx broadening in initial state 
interactions. Although questions have been raised 
concerning issues such as decrease in entropy of the 
system and the spatial extent of the recombining 
subsystems, it is useful to compare the experimen- 
tal measurements with this model to investigate the 
relative importance of various physical processes. 

Fig. |H1 compares model predictions with the mea- 
sured YAsym for TT^ +TT~ +P+P for the rapidity region 
0.5 < \y\ < 1.0 in 0-20% central d-FAu collisions. 
The model predictions from Ref. ^ are consistent 
with the data. Since the pions are the dominant 
hadrons produced in d-|-Au collisions, the YAsym for 
7r++7r^+p-|-p is dominated by them. In the absence 
of predictions from this model for identified hadrons 
in c?+Au data, it is not clear if it can describe the 
YAsym for P+P- One of the reasons for the success 
of the recombination model in the intermediate pt 
range is that in this model a baryon is formed by 
recombination of three shower and thermal partons, 
while a meson needs only two, resulting in a higher 
yield at larger momentum for baryons 0. A com- 
parison of model calculations separately for charged 
pions and p+p is of interest to see if the observed 



E. Comparison to the saturation model 

Finally, we compare our charged pion measure- 
ments to calculations from saturation models 4] . In 
such models the particle production is determined 
by the high gluon density in the Au nucleus and 
the deuteron. The model had successfully described 
the suppression of high px hadron yields at forward 
rapidities for d-|-Au data relative to p+p data at 
RHIC. In contrast to a naive multiple scattering 
picture, where one expects enhancement due to the 
Cronin effect to be more significant for larger for- 
ward rapidities due to the increase in the number of 
scattering centers while probing smaller values of x, 
the saturation models give a completely opposite re- 
sult 28] . For this model the momentum range where 
YAsym > 1 is determined by the saturation and geo- 
metrical scales in the model, as well as the onset of 
the gluon saturation. 

In Fig. El we compare the YAsym data for charged 
pions with the YAsym predictions for inclusive 
hadrons. Such a comparison is reasonable as tt^+tt" 
are the dominant hadrons produced in d-l-Au col- 
lisions. Further, the YAsym values for Tr^-j-Tr" are 
similar to those for p+p. The model calculations are 
in reasonable agreement for \y\ < 0.5 and give the 
correct decreasing trend for YAsym vs. pt- The pre- 
diction of a strong centrality dependence at midra- 
pidity is not observed . Such models are expected 
to work better at forward rapidities at RHIC. The 
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FIG. 10; (Color online) Variation of nuclear modifica- 
tion factor (-Rcp) for tt'^+-k~ and p+p with rapidity for 
PT > 2.5 GeV/c for d+Au collisions at ^/s^ = 200 GeV. 
Also shown for comparison are the _Rcp values for inclu- 
sive charged hadrons as a function of pseudorapidity |^. 
The errors shown as boxes are the systematic errors. The 
error due to number of binary collisions is ~ 14% and is 
not shown in the figure. 



models give larger asymmetries than data for 0.5 < 

\y\ < 1.0. 

V. NUCLEAR MODIFICATION FACTOR 

The gluon saturation effects are believed to mani- 
fest themselves in terms of suppression of transverse 
distributions below the saturation scale. The on- 
set of gluon saturation and the saturation scale, in 
turn, depend upon the gluon density and the rapid- 
ity of the measured particles. The saturation scale 
at RHIC is expected to be ~ 2 GeV^ and depends 
on the colliding nuclei and rapidity as ~ A^l'^e^y 
0,111123. The value of A hes between 0.2-0.3 
and is obtained from fits to HERA data [s^l- It 
is important to study the variation of the nuclear 
modification factor (i?cp) a-s a function of rapid- 
ity. The Rcpiy) and the yAsym(PT) together can 
provide a more stringent constrain on particle pro- 
duction models. Although the present data do not 
have large rapidity span, we will explore the varia- 
tion of Rcp for identified hadrons from forward to 
backward rapidity. i?cp is defined as 



Rcp — 



central 

{dm/dpTdri/{N^i,,))\p,,ipi,' 



where d'^N/dpTdrj is the differential yield per event 
in d-l-Au collisions for a given centrality class. 

Fig. ^1 shows Rcp for tt^+tt^ and p+p with 
rapidity \y\ < 1.0. There may be a decrease in 
Rcp for tt^+tt^ from backward rapidity (gold-side) 
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FIG. 11: (Color online) Variation of tt^/tt^ and p/p 
with rapidity for 2.5 < pt < 10 GeV/c for peripheral 
(40-100%) and n-tag events for d-|-Au collisions at ^snn 
= 200 GeV. Also shown for comparison are the p/p for 
minimum bias p+p collisions. The boxes are the sys- 
tematic errors. The ratios for n-tag events are shifted 
by 0.05 units in rapidity and those for p+p collisions by 
-0.05 units in rapidity for clarity of presentation. 



to forward rapidity (deuteron-side). Within the 
systematic errors (shown as boxes) the Rcp for 
proton+anti-proton is almost constant within the 
rapidity range studied. Also shown for comparison 
are the Rcp values for inclusive charged hadrons as 
a function of pseudorapidity 5]. The dependences 
are slightly weaker than observed by BRAHMS for 
inclusive charged hadrons in the forward rapidity re- 
gion j8ll |. 



VI. PARTICLE RATIOS 

Figure ^2 shows the tt~ /n^ and p/p ratios for 
2.5 < pt < 10 GeV/c as functions of rapidity for pe- 
ripheral (40-100%) and n-tag events for d+Au colli- 
sions at Y^SNN = 200 GeV. The 7r~/7r+ ratio is unity 
for both n-tag and peripheral events in the negative 
(gold-side) rapidity region. For the positive rapidity 
region, the absolute value of ix^ /ix^ ratio is smaller 
for n-tag events compared to peripheral d-l-Au data; 
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FIG. 12; (Color online) Variation of dou- 
ble ratio (7r"/7r+)„_ta9/(7r~/7r+)4o_ioo% and 
ip/p)n-tag / (p/p)4o-ioo% with rapidity for 2.5 < pt < 10 
GeV/c for d+Au collisions at ^fs^ = 200 GeV. The 
boxes are the systematic errors. 



however considering the systematic errors (boxes), 
they are also consistent with unity. The systematic 
errors do not allow for any strong conclusions re- 
garding the differences, which are expected from the 
valence quark and isospin effects at high for n-tag 
events. The p/p ratios are similar within the system- 
atic errors for the two event classes. The p/p ratios 
are slightly smaller than observed for p+p data. 

In order to cancel out most of the systematic 
errors (listed in Table IIII|I . we have plotted the 
double ratio of (7r"/7r+)„_tag/(7r~/7r+)4o_ioo% and 
{p/p)n-tag/{P/p) 40-100% in Fig. US The double ra- 
tio clearly shows the difference between the 7r^/7r+ 
ratio in the forward and backward rapidity regions 
when we compare peripheral d-l-Au collisions and n- 
tag events. The difference for p/p is observed also for 
both rapidities. The boxes shown in the Fig. 1121 are 
systematic errors on the double ratio, which were 
calculated by varying: the distance of closest ap- 
proach of the tracks from the vertex (error of ^1%), 
dE/dx cuts (error of ~1%), px cuts (error of ^2%) 
and small change in rapidity range (error of ^^3%). 
The total systematic error on the double ratio is ~ 
4%. 

The particle ratios and the double ratios can be 
used to get an idea of relative fragmentation of d 
and u quarks to protons, as well as M-quarks to 7r+ 
and 7r~ . Assuming that gluons fragment equally to 
tt"*" and tt^ at high and using the measured dou- 
ble ratio of 7r^/7r+ (Fig. I12|) . one can show that the 
variation of the ratio of u-quark fragmenting to vr^ 
to u-quark fragmenting to 7r+ is sensitive to the frac- 
tion of pions (x^^jgj) originating from quark jets: 



1 




q-jet 



FIG. 13: (Color online) Variation of ^^1,1+ '^i^h fraction 
of pions {Xq_jf,t) originating from quark jets (solid line) 
obtained from the measured ratios of tt^/tt^ in periph- 
eral d-l-Au collisions and n-tag events. The dashed lines 
reflects la systematic errors in obtaining the "^^4. . The 
horizontal shaded band reflects the possible x^_j^f range 
obtained from NLO pQCD calculations using AKK frag- 
mentation functions 13211. 



where r^^ is the double ratio of tt /tt^ shown in 
Fig. ini The results are shown in Fig. [l21 The 
dashed lines reflect the Icr uncertainty in the ratio 

"3^^^ due to uncertainty in the measurements of the 
double ratio of 7r~/7r+. The horizontal shaded band 
reflects the x^-jet value for charged pions from NLO 
pQCD calculations using the Albino-Kniehl-Kramer 
(AKK) set of fragmentation functions (FFs) f33|. 
The width of this band reflects the uncertainty as- 



sociated with 



from the NLO pQCD calcu- 



lations. These are obtained by varying the factor- 
ization scales from 0.5px to 2px- Since the NLO 
pQCD calculations with AKK FFs agree reasonably 
well with charged pion measurements at RHIC, from 
the figure we can conclude that "^^4, ^ 0.3 to 0.6. 

Similarly, it can be shown that the ratio of u-quark 
fragmenting to protons, to d-quark fragmenting to 
protons is given as. 



P 



1 



(1 



I p2 



1.3 ±0.42, 



"q-jet 



where r^i is the ratio p/p ^ 0.81, and rp2 is the 
double ratio of p/p shown in Fig. E| The above 
relation is valid for the assumption that all the anti- 
protons at high px are from gluon fragmentation and 
that gluons fragment equally to protons and anti- 
protons. Similar measurements have been carried 
out by the OPAL collaboration for e"*" -|- e~ colli- 
sions as a function of Xp {2\)/ ^/s, largest scaled mo- 
mentum) ^3 ■ They observed the above ratio to be 
0.637 ± 0.173 ± 0.083 for Xp > 0.2 and that it de- 
creases to 0.165 ± 0.421 ± 0.257 for Xp > 0.5. 
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VII. SUMMARY 

We have presented transverse momentum spec- 
tra for identified charged pions, protons and anti- 
protons from d-l-Au colUsions in various centraHty 
classes at ^snn = 200 GeV. The transverse mo- 
mentum spectra are measured in 4 rapidity bins for 
— 1 < y < 1 over the range 2.5 < pt < 10 GeV/c. 
The rapidity, pt, centraUty, and species dependence 
of the rapidity asymmetry YAsym has been stud- 
ied. We have also presented the rapidity dependence 
of the nuclear modification factor and the 7r~/7r+ 
and p/p ratios for the rapidity range \y\ < 1.0 and 
PT > 2.5 GeV/c. 

The lAsym is found to be larger for 0.5 < \y\ < 1.0 
than for |y| < 0.5 in the range 2.5 < pt < 5.0 GeV/c. 
For higher pT the YAsym approach 1 for both charged 
pions and p+p. From these observations we conclude 
that possible sources of nuclear effects in d-l- Au colli- 
sions, such as parton saturation, nuclear shadowing, 
or energy loss in cold nuclear matter, have a strong 
rapidity dependence which vanishes for > 5.5 
GeV/c. The observed YAsym vs. px dependence 
rules out models based on incoherent initial multiple 
partonic scattering and independent fragmentation. 

Comparison to models based on nuclear shadow- 
ing reveals that incorporation of extremes of gluon 
shadowing at low x does not reproduce the measured 
YAsym- This provides an upper limit on the contri- 
bution of nuclear shadowing to the YAsym- Models 
incorporating multiple scattering, dynamical shad- 
owing, and energy loss in cold nuclear matter are 
in reasonable agreement with the data for |j/ | < 
0.5. However, the Y^sym being independent of 
(> 3 GeV/c) is inconsistent with the measurements 
at higher rapidity. Qualitatively, features of mono- 
tonic decrease in YAsym with px and Rqp with y 
are in agreement with color-glass-condensate (CGC) 
type models. However, there is a lack of quantita- 
tive agreement at higher rapidities where this model 
is expected to work better. Further, the absence of 
very strong centrality dependence at midrapidity in 
the data is in contrast to the predictions from CGC 
models. 

The EPOS and recombination models are in best 
quantitative agreement with the data. The actual 
test of the recombination model is only possible 
when the calculations are available for Y^sym for 
identified baryons and mesons. It will be interest- 
ing to see if this model can explain the observed 
weak species dependence and similar px dependence 
of YAsym for 7T^+TT~ and p+p. 

In general, the study of identified hadron YAsym 
as a function of many variables (y, px, centrality 
and particle type) for d+Au collisions has been able 
to provide some definitive insight on mechanisms of 
particle production in d+Au collisions at y^snn = 
200 GeV. The ^AsymlPx) together with Rcp{y) can 



provide a more stringent constrain on particle pro- 
duction models. It may be mentioned that a detailed 
of study of particle yields (px < 3 GeV/c) at midra- 
pidity and forward rapidity in STAR has revealed 
a possible alternative explanation of the pseudora- 
pidity dependence of Rqp from a purely geometrical 
picture. The decrease in Rqp from negative (back- 
ward) to positive (forward) rapidity can be explained 
by considering the inital asymmetry in particle pro- 
duction in d-l-Au collisions compared to the symmet- 
ric p+p collisions [33| - 

The ratios 7r~/7r"'" and p/p have been studied for 
peripheral and n-tag events for d+Au collisions to 
see the possible valence quark effect. For the range 
2.5 < px < 10 GeV/c and the rapidity region on 
the deuteron-side, the ratios for n-tag events are 
smaller than for peripheral events. However, within 
the systematic errors it is difficult to make strong 
conclusions of valence quark effects on particle pro- 
duction at high px- The p/p ratios are observed to 
be systematically lower than corresponding values 
from p+p collisions. 

The double ratio between n-tag events and 40- 
100% peripheral collision events does reveal a clear 
enhancement in tt"*" production relative to 7r~ at for- 
ward rapidity (deuteron-side). No such enhance- 
ment is observed at the backward rapidity (gold- 
side). Using the above ratio measurements we have 
found ~ 0.3 - 0.6 and ^ ~ 1.3 ± 0.42 for 

2.5 < px < 10 GeV/c. 

A future, high statistics run for d+Au collisions 
at RHIC may be able to provide data that will lead 
to a still better insight into valence quark and gluon 
contribution, as well as isospin effects at high px- 
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